CONCENTRATION PHENOMENA FOR NEUTRONIC MULTIGROUP 
DIFFUSION IN RANDOM ENVIRONMENTS 



SCOTT N. ARMSTRONG AND PANAGIOTIS E. SOUGANIDIS 

Abstract. We study the asymptotic behavior of the principal eigenvalue of a weakly cou- 
pled, cooperative linear elliptic system in a stationary ergodic heterogeneous medium. The 
system arises as the so-called multigroup diffusion model for neutron flux in nuclear reac- 
tor cores, the principal eigenvalue determining the criticality of the reactor in a stationary 
state. Such systems have been well-studied in recent years in the periodic setting, and the 
purpose of this work is to obtain results in random media. Our approach connects the 
linear eigenvalue problem to a system of quasilinear viscous Hamilton- Jacobi equations. By 
homogenizing the latter, we characterize the asymptotic behavior of the eigenvalue of the 
linear problem and exhibit some concentration behavior of the eigenf unctions. 



1. Introduction 

We study the behavior, as e — > 0, of the principal eigenvalue and eigenfunction of the 
weakly-coupled, cooperative elliptic system 



'X 



1.1) - e 2 tr (Aa (-, w) £>V a ) + eb a (-, w) • D<pP a + £ c Q/3 (: 

m 



m) 



13=1 

subject to the conditions 

(1.2) <^ > in U and (p% = on dU (ct = l,...,m). 

Here m > 1 is a positive integer and U C H d is a bounded domain. The unknowns are the 
eigenvalue A e = X E (co,U) and the eigenfunctions (<£%{•, u)) 1<a<m . The underlying random 
environment is described by a probability space (fl, iF, P), and the coefficients A a , b a , c a p 
and a a /s are functions on R d x Q which are required to be stationary and ergodic. (Precise 
hypotheses are found in Section [2] below. ) 

The expectation is that large amplitude, high-frequency oscillations persist as e — > 0, and 
the goal is to describe these oscillations. 
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Problem (II. ip has been proposed and extensively studied in periodic media by physicists 



as a simplified model for the neutron flux in nuclear reactor cores, see [15|, [23j, |24|, |25|, |33 |. 
The modeling assumption is that neutrons are moving in the reactor core (the domain U) in 
m distinct energy groups, each group consisting of neutrons with a similar amount of kinetic 
energy. The function tp e a is the steady-state distribution of neutrons in the ath energy group 
inside the core. The matrix A a describes the diffusion of the neutrons in the ath group, 
the vector b a is the drift, c a/ 3 is the total cross section, which represents the interaction 
of neutrons in various energy groups, and a a p models the creation of neutrons by nuclear 
fission. The factors e 2 and e appear in front of the diffusion and drift terms, respectively, 
due to a physical assumption that the order of the diffusion and drift should be the same 
as that of the microscopic lattice. The principal eigenvalue A £ , in particular whether A £ is 
greater or less than 1, determines the criticality of the reactor. Hence characterizing the 
asymptotic behavior of A £ is of particular importance. We remark that while the model is 
typically written in divergence form, if the matrices A a are sufficiently regular, it can be 
recast in the form of (jl.ip . 



A very complete mathematical analysis o f (II .11) - (1 1.2 1) in the case of periodic coefficients 



was performed by Capdeboscq [l3[ (see also [3lT 13. 4. la. B. I3I l6l 13(1 17|). It was shown in [l3| 
that the eigenvalue A £ admits the expansion 

(1.3) A £ = A + e 2 fi + o(e 2 ) as e ->• 0, 
and the eigenfunctions can be factored as 

(1.4) ^ = ^ a (^)exp(-0.|) + o(l)) as 0. 

Here A G R, 9 e R d , and the periodic function ijj a = ip a {i)) are identified via an optimization 
of a periodic (cell) eigenvalue problem, while (/x, u) is the solution of an effective "recentered" 
principal eigenvalue problem in the macroscopic domain U. Observe that the oscillations of 
the coefficients on the microscopic scale e not only induce oscillations in the solution on a 
scale of e, but also produce a large macroscopic effect, namely an exponential drift. 

The random setting is different. As we will see later, we cannot expect (11.31) and (11. 4p to 
hold in full. Moreover, the approach of [13) does not seem to yield itself to the analysis of (11. ip 
in random environments. Instead, we present an alternative approach. The classical Hopf- 
Cole transformation converts (11.11) into a quasilinear (viscous) Hamilton- Jacobi system, and 
we observe that this nonlinear system may be analyzed by the methods recently introduced 



by Lions and Souganidis [28] and developed further by the authors [9|. Our main result 
is the assertion that there exists a deterministic constant A such that, as e — > 0, A £ — > 
A almost surely in u, together with a characterization of A. In fact, A is identified in 
terms of a convex effective Hamiltonian H. Furthermore, we exhibit concentration behavior 
for the eigenfunctions (p e a , showing that, under some additional hypotheses on the random 
environment, we have, as e — > 0, 

— e\og{p £ a — )• 9 ■ x locally uniformly in U and a.s. in uj, 
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where 9 := argmin^d H . We thereby justify, in random environments, the leading term in 
(II. 3p and, under stronger assumptions, in (jl.4p . 

The article is organized as follows. In Section |2] we state the assumptions and some 
preliminary results needed in the sequel. The main results are presented in Section |3j In 
Section H] we prepare the homogenization of the transformed nonlinear system by studying 
an auxiliary "cell" problem. In Section |5l we define H and do most of the work for the proof 
of Theorem (TJ which is given in Section [6j This analysis is applied to the linear system in 
Section where we prove Theorem [2j In Section [HJ we show that H is strictly convex in p 
(and therefore the eigenfunctions concentrate) in uniquely ergodic environments. 

2. Preliminaries 

2.1. Notation. The symbols C and c denote positive constants, which may vary from line 
to line and, unless otherwise indicated, do not depend on u. We work in the (i-dimensional 
Euclidean space R d with d > 1, and we write R + := (0, oo). The set of rational numbers is 
denoted by Q. The set of n-bj-d matrices is denoted by M nxd , and S d C M, dxd is the set of 
d-by-d symmetric matrices. If v, w G R d , then v <8> w G S d is the symmetric tensor product 
which is the matrix with entries \{viWj + VjWi). For y G R d , we denote the Euclidean norm 
of y by \y\, while if M G M, nxd , M t is the transpose of M. If M G M dxd , then tr(M) is the 
trace of M, and we write \M\ : = tr(M'M) 1 / 2 . The identity matrix is I d . If U C R d , then 
|?7| is the Lebesgue measure of £/. Open balls are written B(y,r) := {x G R d : |x — y| < r}, 
and we set S r := B(0, r). The distance between two subsets U, V C R d is denoted by 
dist(C/,F) = inf{|x — y| : x £ U, y €L V}. If 17 C R d is open, then USC(C/), LSC(C/) 
and BUC(?7) are respectively the sets of upper semicontinuous, lower semicontinuous and 
bounded and uniformly continuous functions U — > R. If / : U — > R is integrable, then we 
use the notation 

l fdv= w\L fdv - 

If / : U — > R is measurable, then we set oscu f '■= ess sup v f — essinf^ /. The Borel cr-field 
on R d is denoted by 23(R d ). If s, t G R, we write s At := min{s, £}. 

We emphasize that, throughout the paper, all differential inequalities involving functions 
not known to be smooth are assumed to be satisfied in the viscosity sense. Finally, we 
abbreviate the phrase almost surely in u> by "a.s. in w." 

2.2. The random medium. The random environment is described by a probability space 
(fl, jF, P), and a particular "medium" is an element wed. We endow the probability 
space with an ergodic group (T y ) y£ - R d of ^-measurable, measure-preserving transformations 
T y : Q — > Q. Here ergodic means that, if D C Q is such that t z (D) = D for every z G R d , 
then either P[-D] = or ~P[D] = 1. An 5F- measurable function / on R d x Q is said to be 
stationary if the law of f(y, •) is independent of y. This is quantified in terms of r by the 
requirement that 

f(y, t z u) = f(y + z, u) for every y, z G R d . 
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Notice that if : Q — > S is a random process, then <f>(y,u) := cj){T y uj) is stationary. Con- 
versely, if / is a stationary function on R d x Q, then f(y, u) = f(0,T y co). 

The expectation of a random variable / with respect to P is written E/, and we denote 
the variance of / by Var(/) := E(/ 2 ) — (E/) 2 . If E G 3, then 1# is the indicator random 
variable for E; i.e., 1_e(cj) = 1 if a; G -E, and = otherwise. 

We rely on the following multiparameter ergodic theorem, a proof of which can be found 
in Becker fiH. 



Proposition 2.1. Suppose that f : R d x — > R is stationary and E|/(0, -)| < oo. JTien 
t/iere is a subset Q C f2 of full probability such that, for each bounded domain V C R d and 
w G fi, 

(2.1) lim / f(y,u)dy = Ef. 

2.3. Assumptions. The following hypotheses are in force throughout this article. The 
coefficients 

A Q : R d x n -> S d , b a : R d x ft -> R d and c a/3 , a a/3 : R d x ft -> R. 
are measurable and we require that 

(2.2) Aq,, 6 a , c a p, and cr a/ 3 are stationary for each l<a,/3<m. 
We assume that, for each u G ft and 1 < a < m, 

(2.3) A a (-,u)eCl;l(R d ;S d ) 
and that A Q has the form A = where, for some C > 0, 

(2.4) E a (3/,w) G M dxm satisfies ||S Q (-, u)\\ c0 ,i (Rd) < C. 

We assume that there exists C > such that, for all u G ft and 1 < a, /3 < m, 

(2.5) H^lctC-, ^)||c , °' 1 (JR d ) + IIM'> UJ )\\c°- 1 (R d ) + lka^(' ; w )llc°^ 1 (R £i ) + ll C a^('' W )llc°^(K d ) - ^ 

The matrices A a are uniformly positive definite in the sense that there exist positive constants 
< A < A such that, for every y, £ G R d , w G ft and index a, 

(2.6) m 2 <A a (y,uj)Z-Z<A\Z\ 2 . 
The matrix is diagonally dominant, i.e., 

m 

(2.7) c a/ 3 < for a 7^ /3, and c Q/ g > 0, 

0=1 

as well as fully coupled in the sense that there exists c > such that 

if {J, 3} is a nontrivial partition of {1, . . . , m}, then for every y G R d and w G ft, 
there exists a G J and /3 <E 3 such that c a p(y,uj) < —c. 
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The hypothesis (12. 8p is satisfied, for example, if c QjQ+ i < — c < for each a G {1, . . . ,m}. 
Finally, we suppose that, for every y G R d , uj G Q, and a, (3 — 1, . . . , m, 

k 

(2.9) a aP > and ^ (T Q7 > C > 0. 

7=1 

We emphasize that (I2.2p -f [2~9l) are assumed to hold throughout this article. 

2.4. Further notation. It is convenient to write the system (11.11) in a more compact form. 
For each a G {l,...,m}, let L a denote the linear elliptic operator which acts on a test 
function ip by 

L a <p := - tr(A a (y, u)D 2 ip>) + b a (y, u) ■ Dip, 
and £ = . . . , L m ) act on $ = (ip u <p m ) by 

£$ := (LxVl, ■ ■ ■ i^m^m) ■ 

The operator corresponding to the microscopic scale of order e is denoted by 

(2.10) (£ £ $)(x):=(£M0(£), 

where ty{x) := <§>{ex). Hence we may write £ e = (L\, . . . , where 

= s 2 tr(A a u?j DV) + ^ (~, w) • ^V- 

In view of the above, the eigenvalue problem (II. ip is written concisely as 

L £ $ £ + Q £ $ e = A e E £ $ £ in U, 
$ £ = on dU, 



(2.11) 



where Q £ = Q £ (x, u) denotes the matrix with entries c aj g(^, u) and E e the matrix with entries 
a a p(~,u)). For future reference, we also write Q = Q\ and £ = Ex. 

2.5. Preliminary facts concerning principle eigenvalues. The full coupling assumption 
(12. 8p endows the linear operators £ and L £ with certain positivity properties related to the 
maximum principle (c.f. Sweers j32j). The Krein-Rutman theorem may therefore be invoked 
to yield the existence of a principal eigenvalue of X £ > of (12. lip , which is simple (has a 
one-dimensional eigenspace) and corresponds to an eigenfunctions $ £ with entries ip £ a which 
can be chosen to be positive in U. We summarize these facts in the following proposition, a 
proof of which can be found for example in (29| . 



Proposition 2.2. Under assumptions (|23j) , fl23D . (1231) . (|27)» , (P and (j2T9]) . i/*e sysiem 

+ Q$ = AiE$ in U, 
$ = on <9£/, 



(2.12) 



has a unique eigenvalue Ai = Xi(U, u) > corresponding to an eigenfunction $ = <&(x,uj) 
with positive entries. The eigenvalue Ai is simple, i.e., the eigenfunction $ is unique up to 
multiplication by a nonzero constant. 
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It is well known (see j29[) that the principal eigenvalue Ai is characterized by the variational 
formula 

(2.13) Xi(U,u) = sup {A G R : there exists * G C 2 (U) k with positive entries 

such that £* + Q* > AS* in U) , 

where the differential inequality is meant to hold entry-by-entry in the classical sense. It 
then follows that the eigenvalue Ai is monotone with respect to the domain, i.e., for each 
oj G Q, 

(2.14) Ai(C/, to) < Ai(V, u) provided that V CU. 

Recalling the scaling relation f 1 2 . X j) between L and £ e , we see that the existence and prop- 
erties of the principal eigenvalue A £ for the problem (12. lip follow from Proposition 12.21 and, 
in fact, 

(2.15) e 2 X £ (U, uj) = Ai (e~ l U, u) . 

Notice from (12.141) and (I2.15P that, if V is any domain which is star-shaped with respect to 
the origin (a property which implies that sV C tV if < s < t), then 

(2.16) the map e \ 6 (V, to) := e 2 \ £ (V, to) is increasing. 

This monotonicity property, combined with the ergodic theorem, implies (see Proposition 17. II 
below) that A e converges, almost surely, to a deterministic constant Ao which is independent 
of the domain U. 

3. Main results 

In this section we formulate our main results, Theorems [T||3] below. To properly motivate 
them, we recall what is known in the periodic setting. To obtain the asymptotics ( 11. 3p and 



( II. 4p . Capdeboscq [12|, |13J introduced the 9 -exponential periodic cell problem 

{m m 
- tr{A a (y)D 2 ij e a ) + b a (y) ■ Di> d a + £ c^(y)^ = \(6) ^ a a p(y)^ in R d , 
j3=l /3=1 
ip a > in R d , y i-)- exp(6 l ■ y)w a {y) is periodic, 

for a parameter 9 G R d . He proved that, as a function of 9, the map 9 h> X(9) is strictly 
concave and X(9) — > — oo as |^| — > oo. This implies that A attains its maximum A at a unique 
9 eR d . Writing 

(3.2) ^( X ) := ^|M an d ^ : = *LlA , 

1 ^ aW <(f) ^ ^ 2 

it was then observed that ( II. ip can be rewritten in the form 

- div(Z(f )Du e a ) = /i e a(f )u £ in U, 



CONCENTRATION PHENOMENA FOR NEUTRONIC MULTIGROUP DIFFUSION 



7 



where the coefficients A and a are periodic and depend on the solution of the ^-exponential 
cell problem (and that of its adjoint). The homogenization of the latter is classical, leading 
to the expansion (I1.3P and factorization (11.41) . 

In the general stationary ergodic setting, we cannot expect a factorization of the form 
(13. 2p to hold in any suitable sense. This is related to the fact that, in random environments, 
correctors do not in general exist (see Lions and Souganidis j26|), and so there is no suitable 
analogue of the functions ip e a . This is not merely a technical problem, and goes to the heart 
of difficult issues in the random setting typically referred to as "a lack of compactness." 

Our analysis in the random case follows a different approach. We begin by introducing 
the classical Hopf-Cole change of variables 

(3.3) ^ e a (x,u) := -elog(p e a (x,u), 
which transforms (11.11) into the nonlinear system 

(3.4) - e tr (A a (f , co) D 2 tp e a ) + A a (f , u>) • D^ a + b a (f , u) ■ 

m 

- ( Cq ^ (f ' u ) - M^K/J (f > w )) ex P ( e_1 « ~ =0 in U. 

(3=1 

The study of the behavior of (13 ,4p as e — > falls within the general framework of random 
homogenization of viscous Hamilton- Jacobi equations. The latter has been studied in the 
scalar case by Lions and Souganidis (13, 2Sl|, Kosygina, Rezakhanlou and Varadhan 20 



and recently by the authors [9(. By adapting the methods of [28| and |9|, we prove a 
homogenization result for (13. 4p . This allows us to understand some aspects of the behavior 
as e of (II. ip and to prove the main result on the asymptotics for A £ and (pf, which is 
Theorem [2] below. 

To explain how the effective Hamiltonian arises, we temporarily "forget" that the eigen- 
value X £ (u) is an unknown in (13 ,4p . This will be accounted for later with the help of 
Proposition 17.11 below. Therefore we consider the system 

^,...,^,/i £ ,f,u;J =g inU, 

where \i e = fi e (uj) > is a (possibly random) parameter, g E C{U) is given, and we define 
H a (p,y,u) := A a (y,u)p -p + b a (y,u) ■ p, 



(3.6) 



f a (zi, . . . , z m , fi, y, u) := ^ (j^afiiv, w) - Cafi(y, u)) exp(z a - z p ) 

13=1 



In writing (I3.5p we have essentially put (13.41) into a more convenient form, replaced A e with 
fi e , and introduced a function g on the right side. 
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Observe that H a = H a (p,y,u) is convex as well as coercive (it grows quadratically) in p, 
while for all £ G R, 

(3.7) f a (z u ...,z m ,p, y, w) = f a (zi + £,-.., z m + f , //, j/, w). 

In addition, there exists C > depending only on the constant in (I2.5P such that, for each 
a — 1, . . . , m, Zi, . . . , z m G R, y G R d , u; G f2 and for all /ii, /i 2 > and pi,p 2 G R d , 

(3.8) |/ a (zi, . . . ,2; m ,/ii,?/,w) - / (*i, . . .,2 m ,jt 2 ,J/,w)| 

< C( max exp(2; a — 2^) ) — /i 2 | • 

and 

(3.9) |i? Q (pi,?/,w) - H a (p 2 ,y,u))\ < C(l + |pi| + \p 2 \)\pi — P2I- 

Our first result is a homogenization assertion for the system (13.51) . We stress that the each 
of assumptions stated in Section [2] is in force in Theorems [TJ [2] and [3j 

Theorem 1. Lei £7 C R d &e any domain, fi £ = yU e (w) a bounded nonnegative random 
variable, and assume that for each e > 0, u G Q and a = l...,m, u e a = u £ a (-,u) is a 
solution of (13.51) . Assume also that there exists u G C(U) and fi > such that, for every 
a = 1 . . . , m, as e — > and a.s. in u, n e — > and u e a — >■ u locally uniformly in U . Then u 
is a solution of the scalar equation 

(3.10) H(Du,p,)=g inU, 

with the effective Hamiltonian H : R d x R + — > R characterized in Proposition I5.il below. 

We will see in Proposition 15.21 that p 1— > H(p, fi) is convex and coercive for each [i > 0, 
while /i H- H(p,fi) is strictly increasing and if(0, 0) < 0. Therefore, we may define A to be 
the largest value of fi for which the graph of p 1— > H(-,fi) touches zero, i.e., 

(3.11) A := sup < > : min H(p, //) < > . 

L pGR d J 

Since if is continuous, 

(3.12) = min H(p,~X). 

We know that, in the random environment, H may have a "flat spot" at its minimum. That 
is, the set 

(3.13) 6 := argmin^(-, A) = {p G R d : ~H(p, A) = 0} 

may have a nonempty interior. However, in certain cases, for example, if p 1 — >- H(p, 11) is 
strictly convex, then p 1— > H(p,X) necessarily has a unique minimum, that is, O = {0}. In 
the latter situation, we obtain that the eigenfunctions exhibit concentration behavior in the 
sense of (13.161) below. 
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We emphasize that A and 9 are deterministic quantities which are independent of the 
domain U. 

We now state the result regarding the asymptotics of the linear system fll.ll) . In what 
follows, U is taken to be a smooth, bounded domain and X e and ip e a together solve the 
system f ll.ip - fll.2p . subject to the normalization 

(3.14) <Pi( x o) — 1 f° r some distinguished x G U and a.s. in uj. 

The following result characterizes the limit of the eigenvalues A e , and uncovers the concen- 
tration behavior of the eigenfunctions ip e a in the case that H{-, A) achieves its minimum at 
a unique point 9. 

Theorem 2. The principle eigenvalue X e (u,U) of (1 1 . 1 j) - ( 1X721) satisfies 

(3.15) e 2 \ £ (u, U) — > A as e — > and a.s. in u, 

where A is given by ( 13. lip . Suppose in addition that if (•, A) attains its minimum at a unique 
point 9 G R d . Then, for each a — 1, ... ,m and as e — > 0, 

(3.16) — £log<^(x, u) — > 9 ■ (x — Xq) locally uniformly in U and a.s. inu. 

We present a sufficient condition for the strict convexity of H. The following theorem 
states that the effective Hamiltonian is strictly convex in p under the additional assumption 
that the random environment is uniquely ergodic. Roughly speaking, this means that the 
limit (12. ip in the ergodic theorem is uniformly with respect to the translations. The precise 
definition follows. 

Definition 3.1. The action of the group {T y ) y& ^d on the environment (Q, P) is uniquely 
ergodic if, for every ^-measurable / : O — > R such that E|/| < oo, there exists a subset 
Q C Q of full probability such that, for every u> G Q, 

i f(T y tu)dy-Ef 

JB(z,R) 

Equivalently, the action of (T y ) y& - R d on the environment is uniquely ergodic if and only if 
P is the unique ^-measurable probability measure which is invariant under the action. 

The space of almost periodic functions may be embedded into the stationary ergodic 
setting (c.f. [8|), and it is easy to see that the resulting random environment must be 
uniquely ergodic. Therefore f)3.16p holds in particular in the almost periodic setting. The 
inclusions are proper: it is well-known that there exist stationary ergodic environments which 
are not uniquely ergodic (for example, any iid environment cannot be uniquely ergodic), and 
uniquely ergodic environments which are not equivalent to translations of an almost periodic 
function [34] . 

Theorem 3. Assume that the action of (r^)^ on (Q, S 7 , P) is uniquely ergodic. Then, for 
each n > 0, 

(3.18) p i y H(p,fi) is strictly convex. 



(3.17) lim sup 



0. 
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Hence = {6}, for some 9 E R d , and the concentration phenomenon (I3.16P holds. 

4. The auxiliary macroscopic problem 

For each fixed S > 0, \i > and p E R d , we introduce the auxiliary macroscopic system 

(4.1) 5v 5 a - ti(A a (y, u)D 2 v 5 a ) + H a (p + Dv 5 a , y, u) + f a (v(, . . . , v 5 m , //, y, u) = in R d , 

with i7 a and f a defined in (13.61) . In the periodic setting, (14. ip is known as the "cell problem," 
and it is central to the homogenization of Hamilton- Jacobi equations. In this section we 
establish the well-posedness of (14.11) . a necessary precursor to the next section, where we 
construct H via a limit procedure using v s . 

Following the usual viscosity theoretic approach, we first give a comparison principle for 
[ITJ. The structural assumptions in Section [2] yield the following proposition, which is 



essentially due to Ishii and Koike |l_9j. The result in [19[ was stated only for equations in 
bounded domains, but we extend it to R d via a simple argument using the convexity of H a . 

Proposition 4.1. Fix 5 > 0, fi > 0, p E R d and tu E Q. Suppose that u a E USC(R d ) is 
bounded above and v a E LSC(R d ) is bounded below, and u = u a and v = v a satisfy, for each 
a = 1, . . . , m, 

5u a - tr(A a (y, u)D 2 u a ) + H a (p + Du a , y, u) + f a (ui, . . . ,u m , fi,y,u) < in R d , 

and 

5v a - tr(A a (y, u)D 2 v a ) + H a (p + Dv a , y, u) + f a {vi, . . . ,v m , fi,y,u) > in R d . 
Then, for every a = 1, . . . , m, u a <v a in R d . 

Proof. According to the structural assumptions, we may select k > sufficiently large (de- 
pending on 5) so that, for each a = 1, . . . ,m, the function (f(y) := k — (1 + lz/1 2 ) 1 ^ 2 is a 
smooth solution of 

dip - tr(A a (y, u)D 2 ip) + H a (p + Dip, y, u) < in R d . 

Modify u a by defining, for each e > 0, 

u a , £ (y) ■= (1 - e)u a (y) + ecp(y). 

Formally, using the convexity of H a and (13.71) . for each a = 1, . . . , m, we have 

5u aj£ - tr(A a (y, u)D 2 u^ £ ) + H a (p + Du aj£ , y, u) + f a (u lj£ , . . . , u m>£ , fj,, y, u) < in R d . 

This can be made rigorous either by using the fact that ip is smooth or by applying 0, 
Lemma A.l]. Since v a is bounded below and u at£ (y) —> — oo as \y\ —> oo, for all R > 
sufficiently large and for each a = 1, . . . , m, we have 

u a , £ < v a in R^ \ B R . 

It then follows from [19|, Theorem 4.7] that, for each a = 1, . . . , m, 

u a ,e < v a in R d , 
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and, after sending e — > 0, the conclusion. □ 

The unique solvability of (14. X)) follows easily from Proposition 14. II and the Perron method. 

Proposition 4.2. For each 5 > 0, fi > 0, p G R d , and u G Q, there exists a unique bounded 
viscosity solution v s (-,ui;p, /i) = (vf(-, p, fj), . . . , ^(-, ui; p, /i)) G C(R d ) m o/ ( 14 -XI) . which 
is stationary. Moreover, there exist C,c> 0, depending only on the constants in ( I2.6p . (12. 5p 
and ( 12. 9p . swc/i t/iat, /or eac/i a; G f2 and a — 1, . . . , m, 

(4.2) - (AH 2 + C(\p\ + fi)) < 5v s a (- 7 u;p,Li) < -(\\p\ 2 - C(\p\ + 1)) in R d . 
Proof. We suppress dependence on u, since it plays no role in the proof. Denote 

:= -i(A|p| 2 + C(|p| + /i)) and ^ a (y) := -^(A|p| 2 - C(|p| + 1)). 

It is easy to check, using ( 12 .6[) . (1 2 . 5 j) . ( 12. 7p and (12. 9j) . that and Tp a are, respectively, a 
subsolution and supersolution of (14.11) in R d . Define, for each a = 1, . . . , m, 

f^(u) := sup {(p a (y) : (fx, . . . , (p m G USC(R d ) are bounded above and 

(p = (if i, . . . , if m ) is a subsolution of ( 14.11) in R d }. 

It is clear from the definition above and Proposition 14 . 1 1 that f a <v & a < Tp a in R d , which gives 
( I4.2p . Standard arguments from the theory of viscosity solutions, utilizing Proposition 14.11 
imply that v 5 G C(R d ) m and v s is a solution of (14. ip . We refer to [14j and to Section 3 of 
[l9j for details. According to Proposition 14. X ^ v s is the unique bounded solution of (14. ip . 
The stationarity of v 5 is an immediate consequence of the stationarity of the coefficients and 
the uniqueness of v s . □ 

In the following proposition, we use the Harnack inequality for linear elliptic systems (see 
Busca and Sirakov jXXj] ) to obtain an estimate, independently of 5, on the difference between 
v s a and vl. The Bernstein method then yields uniform Lipschitz bounds on v 6 . 

Proposition 4.3. For each 5 > 0, > 0, p G R d and u G £1, the unique bounded solution 
v s (-, lj; p, /i) of ( 14. ip belongs to C 2 (R d ) m and there exists C > 0, which depends on upper 
bounds for \p\ and /x but is independent of 5 and u, such that 

(4.3) max ess sup \v^(-, uS) — v t(-, uS) I < C 

Q,/3G{l,...,m} Rd 1 

and 

(4.4) max esssup |.Df^(-,u;)| < C. 

oG{l,...,m} j^d 

Proof. For convenience, we omit the explicit dependence on u since it has no role in the 
argument. The smoothness of v s is immediate from classical elliptic regularity. The estimate 
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(14. 3 p is a consequence of a Harnack inequality for a linear cooperative systems. To see this, 
we observe that w 5 = w s a with w s a := exp(— v s a ) is a classical solution of 

- tx(A a (y)D 2 w s a ) + (2A a (y)p + b a (y)) ■ Dw s a 

m 

- (A a (y)p ■ p + b a (y) ■ p + 5v s a )w s a + ^(c Q/3 (?/) - /la^w^ = in R d , 

0=1 

which is a cooperative, fully coupled linear system, and thanks to (14.21) . has bounded coeffi- 
cients. The Harnack inequality found in [111, Corollary 8.1] yields that 

w 5 a (y) < C w 5 p(y) for each y e R d and a, (3 = 1, . . . , m. 

Rewriting this inequality in terms of v s a and vl yields (14. 3p . 

According to (14.31) . the last term fa(vf , . . . , v^, fi, y) in (14. ip is bounded independently of 
5. We next use the Bernstein method to obtain the Lipschitz estimate (I4.4p . Although it 
proceeds very similarly as the proof of [27|, Proposition 6.11] for the case of a scalar equation, 
for the convenience of the reader we give complete details here because the form of the system 
(14. ip complicates the argument somewhat. For ease of notation we do not display the explicit 
dependence of v 5 a on 5. Select a cutoff function ip e C°°(R d ) such that 

(4.5) 0<<p<l, y = lonSi, if = in R d \B 2 , \D 2 y\ < C(p* and \Dip\ < C<p*. 

It suffices to choose for example (p = ip 4 , where ^ is a cutoff function satisfying the first three 
conditions of (14. 5p . Denote £ Q := \Dv a \ 2 and w a := <p£ a = ip\Dv a \ 2 . An easy computation 
yields 

(4.6) Dv a = £ a Dip + 2ipD 2 v a Dv a , 

(4.7) D 2 w a = £ a D 2 (p + 2Df g> (D 2 v a Dv a ) + <p(D 3 v a Dv a + D 2 v a D 2 v a ). 

Differentiating ( 14. ip with respect to y iy multiplying the result by fv am and using ( 14. 6 p and 
( 14 .7p . we obtain after some calculation that, on the support of if, 

(4.8) 5w a - tr (A a D 2 w a ) + <p tr(D 2 v a A a D 2 v a ) - tpDv a ■ tr(D y A a D 2 v a ) 

+ £ a tr(A a D 2 ip) + A a ^~ l D V ■ {Dw a - 
+ ipDv a ■ ((D y A a (p + Dv a ) + D y b a ) ■ (p + Dv a )) + A a (Dw a - £ a D<p) ■ (p + Dv a ) 

m 

+ \b a ■ (Dw a - £ a D<p) + ipDv a ■ e Va - v W y {na a p - c aP ) 

m 

+ vY. eVa ~ Vf) &Vaf> - c aP ) (\Dv a \ 2 - Dv a ■ Dvp) = 0. 

13=1 
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Now suppose that x G B 2 and a' = 1, . . . ,m are such that 

(4.9) w a /(x) = max sup w a (x). 

ae{l,...,m} xeK d 

To simplify the notation we assume that a' = 1. Then Dw\{x) = and D 2 w\{x) < 0. Using 
these together with (14.3j) . (14.51) and the observation that (14.91) implies that, at x = x, 

m 

(jjutw - c w ) (\D Vl \ 2 - Dvi ■ Dvp) > 0, 

0=1 

after some work we obtain, from (14.81) . 

V ti{AM 2 ) < C(<p\q\(\M\ + 1) + |?| 2 (|L>Vl + 1^1 + ^Dri 2 ) + \q\ 3 (p + \Dp\)) , 

where for convenience we have written M := D 2 Vi(x), q := Dv\(x) and A = A\(x). Applying 
some elementary inequalities and using (12 .4p we get 

ptr(AM 2 ) <C v [ip+ \q\ 2 ^ + |g|V*) + w\ M ?, 

where r\ > is selected below. By (12. 6j) and the Cauchy-Schwarz inequality, 

\M\ 2 < C(tr(AM)) 2 < Ctr{AM 2 ). 
Hence by making rj > small we get 

(4.10) ipti(AM 2 ) < c(l + \q\ 2 <p* + |g|V*) . 
Using the PDE (J4TT]) and the estimates (j4T2l) . we have 

(4.11) (tr(AM)) 2 = (ffi(p + Dvi, y) + /^Vj, . . . , v k , fM, y) + 5 Vl f > X(\q\ 2 - C) 2 . 
Putting (I4.10P and (14.1 II) together, we obtain that 

|g|V < C(l + \q\ 2 ^ + • 
This yields an upper bound on \q\ 4 <f and hence 

( Wl (x)) 2 = |g|V < |5|V < C- 



Thus 



max sup I -D^q, I ~ < max sup I w a I = u^i (£) < C. 

ae{l,...,m} ^ Qg{l,...,m} j^d 



Since the constant C > in the last inequality is independent of the fact we centered our 
ball at the origin, the proof is complete. □ 

Remark 4.4. The essential boundedness and stationarity of v a and (I4.4p imply that 

(4.12) E[£>t£(0,-)] =0. 

This follows from an argument of Kozlov (2lj], see also (l, Lemma A. 5]. 
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We conclude this section by studying the dependence of the solution of (14. ip on the 
parameters p and /x. It is a necessary ingredient in the proof of Proposition 15. II and will yield 
important properties of H. 

Proposition 4.5. Let v s = v^(-,co;p, /i) be as in Proposition \4 ■ 4 Then: (i) for each k > 
there exists C > 0, depending on m, k and the constant in (12.51) . such that, for all 5 > 0, 
uj G ft, Pi,p 2 G H d and < fi < k, 

(4.13) max sup 5\v 5 a (-,u;pi, fx) - v 5 a (-,uj;p 2 , fjt)\ < C(l + IpA + \p 2 \) \pi - p 2 \. 

a£{l,...,m} jjd 

(ii) for each k > there exist C,c> 0, depending only on m and the constants in (12. 9p and 
(12. 5p . such that, for all 5 > 0, uj G VL, p G -8(0, k), < /ii < /i 2 < A; and a = 1, . . . , m, 

(4.14) c{yL2-yLi)<8v 5 a {-,oj-,p,ii 1 )-5v 5 a {-,uj-,p,yL 2 )<C{ii2-iii) in R d . 

Proof. The proof of (I4.13P closely follows the proof of [9j, Lemma 4.6]. We fix pi,p 2 G R d , 
/i > and u; G Q and write v 5 ai {y) := w; Pi, /i) for i G {1, 2} and a = 1, . . . , m. Define 
A := (1 + \pi\ + \p 2 \Y l \p\ - Pi\ and set 

«£(y) := (1 - A)< 2 (y) = (1 - A) ((p 2 - p x ) • y + < 2 (y)) + A(A^(1 - A)(p a - p 2 ) • y). 
It is easy to check, using the convexity of H a and (13. 7p . that w 5 a satisfies 

8w s a - tr(A a (y 7 uj)D 2 w 5 a ) + H a (p + Dw 5 a , y, uj) + f a (w{, ...,w s m ,n, y, uj) 

<XH a (\- 1 (p 1 -(l-X)p 2 ),y,u J ) inR d . 

Since 

A _1 |pi - (1- A)p 2 | < 1 + \ Pl + 2\p 2 \, 
there exists a constant C > depending only on the constant in (I2.5p . so that 

XH a (X^( Pl - (1 - X)p 2 ),y,u) < C(l + | Pl | + |p 2 |)| Pl -p 2 \. 

By subtracting 5 _1 C(l + \p\\ + \p 2 \) \p\ — p 2 \ from u>^, we obtain a subsolution of (14.11) . and 
Proposition 14.11 yields 

(1 - A)< 2 = w s a < < x + 5~ l C(l + M + |p 2 |) | Pl -p 2 |. 

Using (14.21) and rearranging, we obtain 

<a - <i < + N + |p 2 |) |pi - p 2 | 

where C depends additionally on an upper bound for /i. Multiplying by <5 and repeating the 
argument with the indices reversed yields (14.131) . 

The first inequality in (I4.14p easily follows from (14.31) . (12.91) . and Proposition 14.11 Indeed, 
the function v 5 {y, uj; p, /ii) — 5 _1 c(/j 2 — /ii) is a supersolution of (14. ip for ji = /i 2 , for a 
c > with appropriate dependencies. The second inequality follows similarly, since the 
function v s (y, uj; p, jii) — 5~ l C(n 2 — /ii) is a subsolution the same equation for sufficiently 
large C > 0. □ 
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5. Construction of the effective Hamiltonian H(p,fi) 

In the next proposition, which is the analogue of [9, Proposition 5.1], we identify H as 
a limit of 5v s as 5 — > 0, and construct a subcorrector. The argument is based on ideas 
introduced in 28| . 

Proposition 5.1. There exist a continuous function H : R d x R + — > R and a subset Q\ C Q 
of full probability, such that, for every p G R d , /i > 0, R > 0, and a = 1, . . . , m, 

(5.1) lim sup \Sv s a (y, -;p,n) + H(p,fi)\ = inL 1 ^,?), 
5->o y6BH/{ 

and, /or every u G fii, 

(5.2) if(p, /i) = — liminf ^^(O, cu; p, /i) a.s. ma;. 

<5->0 

Moreover, there exists w = w a (y, cu; p, /i) such that, for every (p, /z, a) G R d x R x {1, . . . , m} 
and u G fii, w a (-,u;p, fi) G C 0,1 (R d ) as we// as 

(5.3) Dw a (-, -;p, /i) is stationary and E[Dw a (0, -;p, = 0, 

(5.4) lim |y| io a (y, w) = 

|y|->oo 

and 

(5.5) -tr(v4 Q (w,a;)L' 2 w Q ) + H a (p + Dw a ,y,u) + f a (w 1 , . . . ,w m , fi,y,u) < H(p,n) in R d . 

Proof. Proposition 14.51 allows us to prove the claim for fixed p G R d and /i > 0, and then 
to intersect relevant subsets of Q for a countable dense subset of (p, /i) in R d x R + . We 
therefore fix p and // and omit the dependence on these variables for ease of notation. 
For a = 1, . . . , m, define 

(5.6) vi(y,u):=v s a (y,u)-v 5 1 (0,u). 

The estimates (14. 2p . (14. 3|) and (I4.4p together with the stationarity of v s are sufficient for the 
extraction of a subsequence 5j — » such that, for every R > 0, as j — > oo, 

- 6jV% ->> # weakly-* in L 00 ^ x fi), 
^ w a weakly-* in L°°(B R x fi), 
D53^ ->> -Dw a weakly-* in L°°(B R x fi), 

for a deterministic constant if = H(p,fi) and functions u> Q (-,w) G C 0,1 (R d ). Standard 
arguments from the theory of viscosity solutions (using in particular the convex structure 
of (14. ip and the equivalence of distributional and viscosity solutions for linear inequalities, 
c.f. Ishii |l7|]) yield that, for each a — 1, . . . , m, w — w a is a solution, a.s. in u, of the 
system (15. 5p . We emphasize that in deriving (15. 5p we rely crucially on the convexity of 
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H a (p, y, u>) in p and of f a {zi, ■ ■ ■ , Zk, fi, y, w) in the differences z a — Zp. According to (I4.4p . 
the gradients Dw a satisfy 

(5.7) sup ess sup |-Dw a | < C, 

«e{i,...,m} R d xn 

and they inherit the stationarity property from the sequence Va ■ From (14.121) we deduce 
that 

E\Dw a ] = lim ElDv 8 *] = lim E\Dv 5 J] = 0. 

j— >oo L J j— s>oo L J 

The ergodic theorem (c.f. Kozlov [21J or [§, Lemma A. 5]) yields (I5.4p . 

Having identified H, we must show that it characterizes the full limit of — 8v s a . The key 
step is to use the comparison principle to show that 

(5.8) — H < liminf 5v^(0, u) and a.s. in uo. 

<5->0 

Denote by Q the subset of Q of full probability consisting u for which (15.41) and ( 15. 5 p hold, 
as well as 

(5.9) swp\Dw a (-,u)\ < C, 

the latter condition holding on a subset of full probability by Fubini's theorem and (15. 7ft . 

Fix oj G Q. Choose 5, r] > 0, let 7 > be a constant to be selected below, and define, for 
each a = 1, . . . , m, 

w 5 a (y, co) := w a (y, u>) - (H + V )5~ 1 - 7 (1 + \y\ 2 ) 1 ' 2 . 

Due to (15 -9p and the fact that |-Dwf,| < |-Dw a | + Cj, we see that w 5 = (wf, . . . , w^) is a 
solution of the system of inequalities 

(5.10) 5w 5 a - ti(A a (y, u)D 2 w 5 a ) + H a {p + Dw 5 a , y, u) + f a (w(, . . . , w 5 m , fi, y, u) 

< 5w a - r] + C7 in R d (a — 1, ... ,m). 

Fix a constant r > to be selected below. Choosing 7 := 77/ (2(7) and applying ( 15 ,4p . we 
may estimate the right side of ( I5.10p for \y\ < r by 

(5.11) 5w a - r] + C7 < 5C V + 5r] 3 r - -77. 

Next we observe that by ( 14.21) . the definition of {yf, and our choice of 7, we have 

(5.12) w 5 a -v s a < V a + CS- 1 -cqr on dB r . 

It follows from (15. lip and (I5.12p that by selecting r = C/5i] for a sufficiently large constant 
C > 0, we obtain, for sufficiently small 5 > 0, 

5w a (y) — 1] + C7 < in £> r and — t> * < on 95 r . 

The comparison principle (c.f. [19?, Theorem 4.7]) yields that, for each a — 1, . . . ,m, 

u£<t;* inR., 
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and, in particular, wf,(0) < f„(0). Multiplying this last inequality by 5 and sending 5 — > 
gives 

-H-Ct]< \immf5v 5 J0 } uj). 

5->0 

Disposing of rj > yields ( 15. 8 p for each a; G Q. 

Since — is the L°°(Q) weak-* limit of SjVa (0, u), the reverse inequality of ( 15. 8ft holds 
and we obtain 

(5.13) — = liminf 5f„(0, u;) a.s. in u. 

<5->0 

Now an elementary lemma from measure theory (c.f. [9], Lemma A. 6]) yields that 

(5.14) Sv s a (0,tu) -> -IT in probability and in L^P). 

We now deduce (|5.1I) from a covering argument and the Lipschitz bound (14.41) (see the last 
step of proof of |9|, Proposition 5.1]). □ 

We next collect some elementary properties of H. 

Proposition 5.2. The effective Hamiltonian H : R d x R + — > R has the following properties: 

(i) for each p G H d , the map fi (->■ H(p,/i) is strictly increasing; 

(ii) for each \i > 0, the map p (->• H(p,fi) is convex; 

(iii) there are positive constants c, C > 0, depending only on the assumptions, such that 

(5.15) X\p\ 2 - C(l + \p\)< Hip, ^)<h\p\ 2 + C(\p\ + n). 

Proof. It is immediate from (14. 14j) and (15. ip that, for all p G R d and < \i\ < yU 2 , 

c(/i 2 - /ii) < #(p, /i 2 ) - H(p, /ii) < C7(/i 2 - /ii) 

for C, c > depending only on upper bounds for \p\ and /i 2 . This yields (i). 

To prove (ii), fix p\,p 2 G R d , > 0, to G f2 , set g := \{pi + p 2 ) and, for each 5 > 0, 

(5-16) := + -v 5 a (y,u;p 2 ^). 

The convexity of f a in the differences z a — zp and the convexity of H a in p easily yield that 
w 5 satisfies 

(5.17) 6w 5 Q - tv(A a (y : co)D 2 w s a ) + H Q (6 + Dw s a , y, 00) + f a « ...,<, n, y, u) < in R d . 

Proposition 14.11 implies that w 5 a < v^(y, u; q, /1). Multiplying by —5 and passing to limits 
with ( 15. ip in mind yields (ii). 

The bounds (I5.15P are immediate from (14. 2p and ( 15. ip . □ 

Remark 5.3. Notice that ( I5.15P implies that H(0, 0) < 0. It follows, then, from (i), (ii) and 
(iii), above, that the number A > given in ( 13. lip is well-defined and 

minIT(-, A) = 0. 
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6. The homogenization of the Hamilton- Jacobi system 

The L l convergence in the limit (15. ip can be upgraded to almost sure convergence. That 
is, we claim that there exists an event Q 2 Q ^ of full probability such that, for every R > 
and ui G fl 2 , 

(6.1) lim sup \5v s a {y,u;p, /x) + H(p, /i) = 0. 

<5^o y eB R/s 

To prove this, the subadditive ergodic theorem must be applied to an appropriately chosen 
subadditive quantity. 

Here we outline a proof of (16. ip which follows closely the ideas of j9[ . Due to the similarity 
to jsj, we omit the details. In fact, the argument is much simpler here since the system is no 
more complicated than the scalar case and, unlike jsj, we are in the context of a bounded 
environment. 

Sketch of the proof of (16. ip . For fixed p G R d and fi > 0, we consider what we call the metric 

problem, which is the system of equations 

(6.2) 

- tv(A a (y, u)D 2 ml) + H a (p + Drri^, y, u) + f a {mj, m^, ft, y, u) = 7 in R d \ B(x, 1), 
coupled with the conditions 

(6.3) x, to; p, fi) = ondB(x,l) and liminf |y| _1 m^(?/, x, cj; p, ji) > 0. 

|y|— >oo 

Here 7 6 R is a parameter, and it is possible to show that (I6.2l) - (l6.3p is well-posed, i.e., 
there exists a unique solution provided that 7 > H(p,fi). In fact, for such 7 there is 
a comparison principle for the system (16. 2 p in exterior domains under very general growth 
conditions at infinity (and see Proposition 6.1 in which is easily generalized to the weakly 
coupled system). 

An argument very similar to the proof of f 14 . 3 1) gives the estimate 

(6.4) max sup x, lo) — mZ(-, x, u) I < C. 

a,/3G{l,...,m} Rd 

The comparison principle then implies that the m^'s are increasing in 7 and jointly stationary 
in the sense that, for every x,y, z e R d and u G f2, 

(6.5) ml(y,x,r z u) = ml(y + z,x + z,u), 

and that, up to a deterministic C > 0, the m^'s are almost subadditive, i.e., for all x,y,z G R d 
and oj G Q, 

(6.6) ml(y, x, u) < ml(z, x, u) + ml(y, z, u) + C. 

The multiparameter subadditive ergodic theorem (c.f. Akcoglu and Krengel [lj) then yields 
that, almost surely in u, 

(6.7) rffiy — x) = lim -mZ{ty,tx,u) 

t— >oo t 
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for a deterministic function WC which, due to (16. 4p . is independent of a. In fact, we can select 
a single event ^ of full probability on which the limit (16. 7p holds for every w G O2, 
p G R d , /i > 0, and 7 > Tl(p, /i). 

With the help of what we have obtained already in Proposition 15.11 we can characterize 
the limit function m 7 . We take a subsequence of #'s along which the convergence in (15. ip 
holds almost surely and argue with a reverse perturbed test function argument (introduced 
in jl, Proposition 6.9]) that 

(6.8) Hip + Dwf,^) = 7 inR rf \{0}. 

Having identified an almost sure limit in terms of the effective Hamiltonian H, we may 
conclude the proof of (16.11) for every u G Q2 by using a perturbed test function argument 
very similar to the one in the proof of [9_|, Proposition 7.1] (or the one below). □ 

With (16. ip in hand, we present the proof of Theorem [U 

Proof of Theorem [IJ We argue only that u is a subsolution of ( 13 . 10[) in U, the verfication 
that it is a supersolution following along similar lines. The proof is by the classical perturbed 
test function method of Evans [161 ]. 

Assume that for some (p G C°°(U) and xo G U, 

(6.9) x h-> (u — tp){x) has a strict local maximum at x = X\. 
We must show that 

(6.10) H{D ( p{x 1 ), fJl )<g{x 1 ). 
Suppose on the contrary that 

(6.11) 7] :— H(D<f(xi), /j,) — g(x 1 ) > 0. 

Fix u G for which fJ, £ (co) — > fi and u £ (-,u) — > u uniformly in a neighborhood of x\. Set 
p = Dip(x\) and define the perturbed function 

where v e a is the solution of (14. ip with S = e. We claim that, for sufficiently small r, e > 0, 

-e tr (A a (l , co)DY a ) + H a (D<p%, f , u) + fj^, . . . , ^, /i e , f , > ^?(x) + ^ in r). 

Indeed, this follows from the continuity of H Q , f a and g, and (13. 7p . (I4.ip . (16. ip and (16. lip . 
The maximum principle for the cooperative system then implies that 

max max (u £ a - <p B ) = max jnax « - cp e a ) . 

ae{l,...,m} dB(xi,r) ae{l,...,m} B(x 1 ,r) 

Using (I4.14p and (16. ip . we send e — > to deduce that 

max max (u a — ip a ) = max jnax (u a — tp a ) . 

aE{l,...,m} dB(x!,r) a€{l,... ,m} B(x!,r) 
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This contradicts (16. 9p for small enough r > 0. We have verified (I6.10p . which confirms that 
u is a viscosity subsolution of (I3.10p . The proof that u is also a supersolution is argued along 
similar lines. □ 

7. Concentration phenomena 

The proof of Theorem [2] is presented in several steps. First, in the next proposition, we 
use Egoroff 's theorem, the ergodic theorem, and the monotonicity in ( 12.161) to show that the 
eigenvalues X £ (U,u) converge almost surely in u to a deterministic limit Ao, which is inde- 
pendent of the domain U. Comparing the eigenf unctions ijj e a and the approximate correctors 
v s a allows us to conclude that A = A, from which the concentration of the eigenf unctions 
follows easily if 9 can be defined unambiguously. 

Proposition 7.1. There exist a subset Qq C Q of full probability and a constant Ao G R 
such that 

limA £ (£7, uj) = Ao for every bounded domain U CH d and all tu G Qq. 

Proof. According to (I2.16p . for each fixed u G Q, the eigenvalue X £ (Bi,u>) is increasing as a 
function of e. Therefore, for every uj G Q, there exists a number Xq(oj) G R such that 

(7.1) X E (B 1 ,u) 4 Xq(u}) ase|0. 
We claim that, for each \i G R, the event 

A r ={wGfl: A (w) > fi} 

has probability P[A M ] G {0, 1}. This follows from the ergodicity assumption once we show 
that t z {Ap) = A M for every z G R d . Indeed, for \z\ < e^ 1 , we have, by stationarity, 

X e (B 1 , r z u) = X 1 (B{z, e- 1 )^) > X 1 (B(0, 2 e -\ u) = X e/2 (Bi,u), 

and similarly, for such e, we also have X £ ^ 2 {Bi 1 t z uo) < X £ (Bi,u). Hence Xq{t z uo) = Xq{uS) for 
every z G R d . This implies t z (A m ) = A M for every z G R d . 
It is then immediate that (17.11) may be improved to 

(7.2) A £ (.Bi, u) I A as e I for every u G f^, 

for some deterministic constant Ao and subset Oi C O of full probability. 

Using Egoroff's theorem we find a subset E C Q with probability P[E] > | such that 

X £ (Bi, uj) l Ao as £ l uniformly in uj G E. 

For each fixed wefl, define the set 

A w := {y G R d : T y uj G E) . 

By the ergodic theorem, for each bounded domain V C R d , there exists Qy Q ^ of full 
probability such that, for each uj G fly, 

lim£l,„Q^=P[£]>i. 
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Choose a countable basis 25 for the Euclidean topology on R d consisting of balls, let f2 2 : = 
flygs^y and define fl Q := Q\ fl Q 2 - 

Fix now a domain U C R d , a small constant 77 > 0, and select an element 6 B with 
V C Z7 and set 5 := dist(V, dU). It follows that, for every u G Qq, there exists T = T (uj) > 
sufficiently large so that, for all ui G E and < e < T (o;) _1 , 

A E (5i, w) - A < r] and eA^nV ^ 0. 

Now fix u G fio. Suppose that < e < 5T (u;) _1 and select y G eA^HV. Then -B(|/, 6) C.U 
and from the stationary hypothesis as well as (I2.14p . (" 12 . 1 5[) and the above properties, we 
may deduce that 

\ e (U, u) < X e (B(y, S),u) = X £/6 (B U truj) <X + r]. 

It follows that 

limsup \ £ (U, oj) < X . 

Owing to the fact that U C B R for some large R > 0, and that u G fi C fi 1; we use (12.1 4p 
and (12. 161) to conclude that 

X e (U, u) > X £ (B R , u) = X £/r (B u uj) > A . □ 

Next we use Theorem [1] to show that Ao equals A defined in ( 13. lip . From this we conclude 
the concentration (13 . 16[) of the eigenfunctions and complete the proof of our main theorem. 

Proof of Theorem^ Let ip £ a denote the functions defined in (13.31) for a = l,...,m and 
normalized according to ipl(xo,uj) = for some fixed xq G U. An argument very similar to 
the one in the proof of Lemma [4.21 yields, for each V CC U, the bound 

sup sup|V£(-,w) < Ce, 

a,/3e{l,...,m} V 

and then the local Lipschitz estimates 

sup sup \Dij) E a {-,u)\ < C, 

ae{l,...,m} V 

for a C > independent of e. Taking a subsequence, also denoted by e, we find ip G C^(U) 
such that, as s y and for every a — 1, ... ,m, 

(7.3) ip e a — > ip locally uniformly in U. 

Now Theorem [1] and Proposition 17.11 imply that ip satisfies the equation 

H(Dip, A ) = in U. 

It follows at once that Aq < A. 
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To obtain the reverse inequality, we select (p,p) such that H(p,fi) < 0. Set 5 > 
sufficiently small so that the event mijj/s Sv^(-, cu; p, /i) > has probability at least |. For uj 
belonging to this event, and, if in addition \ s (U,u) < /i, the map 

x4 min (tp 5 a (x,u) ~ v s a (5x,u;p, //)) 

a£{l,...,m} 

cannot have a local minimum in U according to the comparison principle. This is a con- 
tradiction, since ip £ a (x,uj) — >■ +oo as a; — >• 9t/ and ?;*(•, cu; p, //) is bounded. We deduce that 
P[Xs(U,uj) > fj] > 5 for small 5 > 0. According to Proposition 17.11 we have Ao > fi, and 
hence Ao > A. Therefore Ao = A and we obtain the limit (13.151) . 

Finally, in the case {p : H(p, A) = 0} = {9}, we obtain that Dip = 9 almost everywhere 
in U . It follows that ip(x) — 9 ■ (x — x ) for each x G U, and hence the full sequence ip £ a 
converges to ip. The concentration behavior (I3.16P then follows. □ 

8. Strict convexity of H(p,p) in uniquely ergodic environments 

We prove Theorem [31 Throughout this section, we assume that the action of {T y ) yG ^d on 
the environment (Q, IF, P) is uniquely ergodic (see Definition 13. ip . 

It is worth revisiting the proof of the convexity of p t- >■ H(p,fi) (Proposition 15.2( h)) to 
see if there is extra information we discarded. The argument essentially comes down to the 
derivation of ( 15.171) . There is no doubt that any strict convexity on the part of H must 
be inherited from the H a 's, which satisfy, for every p±,P2 £ P d with q := ^pi + ^p 2 and 

(y, uj) G R d x Q, 

(8.1) -H^p^y^u) + ^H a (p 2 ,y,u) - H(q,y,u) = i(pi - p 2 ) ■ A(y,u)(p l -p 2 ) 

Using ( 18. ip . we observe that, with w 6 a defined as in (I5.16p . we may improve ( I5.17P to 

(8.2) 5w s a - tr{A a {y, co)D 2 w 5 a ) + H a {q + Dw s a , y, u) + f a (w(, ...,w 5 m ,fx, y, u) 

< ~-X\pi -P2 + Dv 5 a (y,u]pi^) - Dv s a (y,u);p 2 ,p,)\ 2 =: -h a (y,u) in R d . 
The hope is to use the term —h a on the right side of (18.21) to show that, for some c > 0, 

wi-viiy^e,^)^^- 1 . 

If h a is bounded below by a positive constant, then the desired conclusion is immediate. 
However, for pi close to p 2 , there is not a definite reason why this should be true. All we 
can say is that h is stationary, nonnegative, bounded, and satisfies, by Jensen's inequality 
and dlllD, 

(8.3) E/i a (0,-)>^Ab 1 -p 2 | 2 . 
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However, what we actually need is something weaker than for h a to be bounded below by a 
positive constant. As we will see, it is enough to rule out the presence of large "bare spots." 
That is, we need to ensure that, for some R > 0, the set on which h a is greater than some 
positive constant takes a uniform proportion of each ball of radius R. This is precisely what 
the unique ergodicity hypothesis gives us, as we see in the following lemma. 

Lemma 8.1. Suppose that g = g(y,u) is stationary, nonnegative, and does not vanish a.s. 
in uj. Then there exist constants r],p > 0, depending on the distribution of g(0, ■), and a 
subset fliCfi of full probability, such that, for every u G fli, there exists R > such that 

(8.4) inf \{y € B(z,R) : g(y) > V }\ > p\B R \. 
Proof. Let 

E:={uen : \{y G £(0,1) : g(y,u) > rj}\ > p|£?i|}, 

with i],p > chosen small enough so that P[E] > 0. According to (13.171) . there exists a 
subset Q C Q of full probability such that, for every u> G £1, there exists R > 1 sufficiently 
large that 

inf / t E {r y co) dy > Jp[£] > 0. 

JB(z,R) 2 

That is, for each u G Q, there exists R > 1, depending on u, such that 

(8.5) \{y G B(z,R) : r y u G E}\ > Cl \B R \, 

with Ci := ^Pf-E 1 ] > 0. Freeze u G fl for the remainder of the argument, let 

(8.6) D(z) := {y G B{z, R) : T y u G E} 
and observe that the stationarity of g yields 

D(z) = {yE B(z,R) : \{x G B(y, 1) : g(x,u) > V }\ > p}- 

According to the Vitali covering lemma, there exist yx,...,y# G D(z) such that the balls 
{B(yi, l)}f = i are disjoint and 

i 

(8.7) ^)C|J%,3). 

i=l 

Since the balls {B(yi, l)}f =1 are disjoint and r y .uj G E, we have, for any z G R d , 

i 

{yeB(z,2R):g(y,u)> V }\ > \ \J{y e B(y t ,l) : g(y,u) > rj}\ > pi^. 

i=l 

It follows from (ESJ and (EZJ) that l\Bi\ > c\D(r)\ > cci\B R \ > c\B 2R \ and hence (lOl . □ 
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In a uniquely ergodic environment, we can prove that many limits derived from the ergodic 
theorem are uniform. Another example is the following useful lemma, which generalizes the 
existence of "approximate correctors" used by Ishii [18j to prove homogenization of Hamilton- 
Jacobi equations in an almost periodic environment (the proof we give below also works in 
the non- viscous setting for first-order Hamilton- Jacobi equations). 

Lemma 8.2. Under the uniquely ergodic assumption, the convergence (16. ip can be improved 
to 

(8.8) 5v^(y,u;p, fi) — >■ H(p, fi) uniformly in R d and a. s. into. 

Proof. Let e > 0. The unique ergodicity assumption and (15.1 ft yield some 8\ > small and 
R > 1 large such that, for all < 5 < 5 , 

sup sup inf \5v^(y, oo; p, p) + H(p, p) | < e. 

a£{l,...,m} y^W B(y,R) 

But then the Lipschitz bound (14.41) gives, for every < 5 < <5oniin{l, (CR)^ 1 }, 

(8.9) sup sup 1 5v s a (y,u;p,p) + H(p,p)\ < s + SCR < 2e. □ 

ae{l,...,m} y eR d 

Next we discuss a so-called "growth lemma," an important analytic tool in the proof of 
Theorem [3j It is a quantitative strong maximum principle which measures how the negative 
term —h a on the right side of ( 18.21) forces w 5 a to be lower in comparison to u; 9, fx). We do 
not give the proof here, since going into details would take us very far off course. However, 
the proof is nearly the same as the proof of the classical growth lemma (c.f. Theorem 2 
on page 118 of Krylov [22]), which follows from the ABP inequality. For this purpose we 
need the following ABP inequality for weakly coupled elliptic systems proved by Busca and 
Sirakov (llj . 

Lemma 8.3. Fix \i > 0, p e H d , uj G Q and suppose that a, r G R and u a and v a satisfy 

(8.10) -tr (A a D 2 u a ) + H a (p+Du a ,y,u) + f a (u 1 , . . . ,u m ,ii,y,u) < r+a-h a (y,u) in B 2R 
and 

(8.11) -ti(A a D 2 v a ) + H a (p + Dv a ,y,u) + f a (v 1 ,...,v m ,fi,y,uj)>T in B 2R , 
where for some r], p > 0, 

(8.12) min |{x G B R : h a (x,oS) > rj}\ > p. 

a£{l,...,m} 

Then there exist constants k, o~q > 0, depending on the constants in the assumptions as well 
as R, 7] and p, such that o < oo implies that 

min inf (v a — u a ) > k + min inf {v a — u a ). 

Q6{l,...,m} Br ae{l,...,m} B2R 

We combine the preceding Lemmata into a proof of Theorem [31 
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Proof of Theorem^ We select /i > 0, p\ ^ p 2 and set q := \{j>i + p 2 ). We argue by 
contradiction under the false assumption that H(pi, p) + H(p 2 , //) = 2H(q, p), proceeding by 
way of a comparison between the functions w 5 a defined in (15.161) and the solutions •; g, p) 
of (14. ip with p = q. 

Fix e > very small and > very large. According to Lemma E21 we may choose 5 > 
sufficiently small to ensure that, for each a = 1, . . . , m, 

(8.13) 8v s a {y,co;q,p)<-H{q,p)+£ in R d 
as well as 

(8.14) w) > -H(q,fj) -e inR d . 
Therefore, we have 

(8.15) - tr {A a D 2 v s a ) + H a (q + Dv s a , y, u) + f a (v(, ...,v s m ,/i, y, u) > H(p, fi) - e 
and, for h a defined in (18.21) . 

(8.16) - tr (A a D 2 v 5 a ) + H a (q + Dv s a , y, u) + f a (v(, ...,<,//, y, u) < H{p, /i) + e - h a (y, to) . 
According to Lemma 18.11 and the growth lemma, if we choose e > sufficiently small, then 

(8.17) M(y):= min (v s a (y) - w 5 a (y)) 

a£{l,...,m} 

satisfies, for some k, R > 0, 

(8.18) M(y)>K+ inf M(z). 

zeB(y,R) 

Such a function cannot be bounded. Indeed, if M were bounded, then, for any > 0, the 
function z i— > M(z) + /3\z\ would achieve its global minimum at some point y e R d . But 
then we would have 

M(y) < M(z) + P(\z\ - \y\) for all z e R d , 

which is incompatible with (I8.18P if we take (3 < k/R, We conclude that M is unbounded. 
However, in light of its definition (I8.17p . the unboundedness of M contradicts the bounded- 
ness of v s a and w s a . This completes the proof. □ 
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